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Abstract. This study presents a generalization of the Landau hydrodynamic solution for multiparticle 
production applied to non-central relativistic heavy ion collisions. The obtained results show longitudinal 
scaling of elliptic flow, V2, as a function of rapidity shifted by beam rapidity {y — ybeam) for different 
energies {^snn = 62.4 GeV and 200 GeV) and for different systems (An- An and Cu-Cu). It is argued that 
the elliptic flow and its longitudinal scaling are due to the initial transverse energy density distribution 



^ and initial longitudinal thickness effect. 
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^ 1 Introduction dinal scaling of elliptic flow [2 . For example, the AMPT 

(a multiphase transport model, ver 1.11) [3] and UrQMD 
Experimentally observed azimuthal asymmetries of parti- (ultrarelativistic quantum molecular dynamics, ver 2.3) 
cle production in non-central heavy ion collisions are cur- g] do not reproduce experimentally observed scaling of 
C3 rently of high interest, as they provide more information elliptic flow. While, the AMPT with string melting and 
about the early dynamics of the high-energy nuclear reac- the Buda-Lund model, based on an analytic solution of 
tions. Moreover, RHIC data [I] of elliptic flow for differ- perfect fluid dynamics, can reproduce experimental data 
ent pseudorapidities shows universal behavior for different |5] ^ 
nuclei and for different beam energies. Not all of the theo- 
retical models are able to reproduce the observed longitu- Landau's approximate hydrodynamic solution for par- 

tide production in relativistic heavy ion reactions [Sl[7] 
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tention and reproduces well the rapidity spectrum of ob- are presented in section |3] The last section presents initial 

served particles in the relativistic nuclear collisions [5]. In state and obtained results. 

recent studies [^llOj the particle production function for The presented solution is approximate, but in compar- 

different rapidities was slightly modified compared to the ison to the computational hydrodynamics is analytic and 

original version to the following form: transparent. The main assumptions of the model coincides 

with original Landau approximations and are as follows: i) 

dN/dyoccxp^iyl^^^-y^), . . 

longitudinal and transverse parts of hydrodynamic equa- 

where ybeam = ln(^SArjv/mAr) is the beam rapidity. In or- tions are solved separately; ii) the equation of state of ideal 

der to make a comparison with experimental particle spec- relativistic gas, P ~ e/3, is used to solve transport equa- 

tra one has to normalize the distribution with the total tions; iii) transverse expansion does not include initial flow 

number of particles, which is unknown from the Landau and is pressure gradient driven. 

solution. The key reason lies in the fact that the Landau 

solution does not conserve total energy, thus the total en- 2 Longitudinal Expansion 

tropy of the system cannot be obtained which leads to an 

In this model we solve the equations of energy-momentum 

unknown number of produced particles Ntot- Anyway, the 

conservation: 
gives the right shape 

9^r^'' = 0, (3) 

of the rapidity distribution for different relativistic ener- 
gies. On the other hand, from the definition of the elliptic where the energy-momentum tensor reads as: 

T^^ = (e -t- P)u^u'' - P.g^"'. (4) 

^^(^^ ^ J d<j>{dN/d^dy) cos(20) 

J d(j>{dN / d(j)dy) solution of local conservation laws (eq. [s]) with the 

one can easily see that normalization of total particle pro- , . r ^ ^ ^ ^^ ^ • r 

t- t- equation 01 state represents the dynamics 01 the system 

duction is not needed. So the task is to modify the approx- , 1 1 n i- r 1 

by relating bulk properties ot the matter, such as: energy 

imate solution of Landau, including the transverse angle , , , „ i.in n r,in-i 

° ° density, e, local pressure, F and the lour-now of the nuid, 

(b dependence in the solution. ,, n/i \ rm r 1 j j -1 

u^^ — u^[l,v). Lhe equations ot the hydrodynamic lon- 

In this study the analytic solution for relativistic hy- j- 1 • • 1 , 1 j- • i 

J J J gitudmal expansion m 1+1 dimension, along the z axis 

drodynamic equations will be presented. In section [2] the j-g^ds 



as: 



dT 



longitudinal expansion will be summarized for the sake of r,^nn .^mor ^mn 

01 01 01 

-lu- + = 0, — — + = 0. (5) 

completeness of the study, even though the same results ci^ 01 az 

can be found elsewhere. The transverse part of the solu- The details on how to solve the above equations can be 

tion with asymmetrical pressure gradient driven expansion found in [6] , [7] , [9] , |10j , |13] , [14] , so here we present the 
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derivation shortly including only main equations and the The transformation back to the (t, z) coordinates is shorter 
result, as we will need it later. with the following new variables: 

Solution of the equations of hydrodynamics ([s]) starts 

y+=\n{{t + z)/A), y^=\n{{t-z)/A), (10) 

by transforming relativistic velocity field components to 

rapidity terms, as: where A is the initial thickness of the system in the beam 

direction, z. Also, A is the initial condition after which 

It" = coshy, = sinhy. (6) 

the equation of state is assumed to be valid and evolution 
From the above transformation, naturally follows, that: equations ^ are applied 

V, = tanhy and (u^-)^ - {u'^f = -1. In order to solve ^^^^ g^^^j ^^^^^-^^ ^^^^^^^ ^^^^-^^^ e{y+,y^), and 

hydrodynamic equations ([s]), the following variables are rapidity y{y+ ?/_) is- 

introduced: 

e-{y+,y~) = eoexp[-4/3(j/+ +y_ - ^y+y^)], (11) 

K = Inw, w~e+p, X = ^ wUzZ — wuqI, 

y(y+,y_) = (y+-y_)/2, (12) 

where x is called the Khalatvikov potential, tp is hydrody- 

, ,. , 1 • 1 • J r. 11 ,1 1 ,• where z = i tanhw. The above solution of 1+1-dimensional 
namic potential, which is denned by the relation: wui — " 

drh 1 / • r c J- i J A J relativistic hydrodynamics equation (|5l will be connected 

where -0 is a function of coordinates and time. And, j j ^ u 

n . ,1 ,1 , TT • 1 i- f J to the solution of transverse expansion, in order to obtain 

w = e + F IS the enthalpy. Using the relation tor sound ^ ' 

, . , multiplicities of produced particles for different rapidities, 

velocity, as: t- t- t t- 

n dw 



c2 = 

s o J 

w an 



where = 1/3 for the ideal relativistic gas equation of 3 Transverse expansion 

state [TS] . After the Legendre transformation to the hodo- 

In order to solve the transverse part of hydrodynamic 

graph plane, the Chaplygin equation for supersonic expan- 

equations ([s]) we will follow original Landau assumptions 

sion reads as: 

with some modifications. For simplicity, polar coordinates 
^ '^sJ ~ ^ ' will be used, where four-flow components in polar coordi- 

, ., , . ,. , ,. r,, 1,-1 nates are: = cix.M, Mo = (1 - (i^J-^ + ^'^i0^))"^/^ 

Detailed investigation on the solution of the 1+1 dimen- "-^ "-^ 

. Ur = uoVr and energy-momentum tensor (HI) components 

sional perfect nuid hydrodynamics can be found m [1611171 ' 

rrm m, , • r , • , ^rc as foUows: 

118) . ihe solution for above equation reads: 

t = e-«r^coshy-|^sinhy), (8) = (e + P)(.") + P, T^^ = (e + P)(.Y.^ + P/.^ 



V Ok dy 



( (^X 9x \ 
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Assuming that the transverse velocity is radial, — 0, the an important role, as it involves transverse asymmetry to 
hydrodynamic equation ([s]) for the transverse dynamics at the solution by making a different transverse pressure gra- 



the fixed transverse angle becomes: dient for different (f> angles. Now from eq. (14), we express 



f)rpOr pirprr the trausvcrsc displacement dependence on time, as: 

Ko^e^*^. (16) 

Inserting energy-momentum tensor expressions in the above v"*^ > 

equation and using ideal gas equation of state, P = e/3. The last stage of the model is the so called conic-flight 

one gets: stage, where energy and entropy fiuxes stop changing for 

fi 2 f) the fixed cone element 2'Krdr. This stage coincides with 

4e("Y^+4e(z.Yf^ + ^=0. (14) 

the kinetic Freeze-Out (FO), as the model does not include 
Following the original Landau derivation, the fist term any hadronic re-interactions after the evolution stops. With 
in the above equation is an acceleration dependence, and the help of the formula ([T6| we obtain a hypersurface in 



is assumed to be equal to dvr/dt = 2r{t)/f. The sec- space-time, after which hydrodynamic fiow stops and mat- 

ond term is set to zero, Vr being comparatively small. To ter streams freely towards detectors. Here again we follow 

simplify the third term. Landau used de/dr sa -c/Ra, the original Landau model, assuming a fixed transverse 

because the energy density at the center has value e and distance, r{tFo) = a = 2i?^, for the conic-flight to start 

is zero at the edge of the system, r ^ Ra- In the case at the distance of the nuclear diameter. Thus, we obtain 

of peripheral collisions, we do not expect centrally sym- the value for the time, when the conic-flight starts, which 

metric energy density distribution, thus the assumption is reads as: 

modified, as: tpo = 2cosh?/J- — (17) 



de _ e(r = R^,) - e(r = 0) 



(^25) where the relation = coshy was used. 



dr i?0 

The solution in the conic-fiight stage is straightfor- 

where R^ is a transverse radius of the system, which changes 

ward, as the energy and entropy does not change at a 

with the angle, as the system is not centrally symmetric. 

fixed cone element. The transverse and longitudinal so- 

We do not know the value of e{R^), so we introduce a new 

lutions are matched at the time t = tpo- Knowing that 

function, f{Rcj,) — e{r — R^)/e{r = 0), which is a fraction 

dS ~ svPdz at a given time within element dz, and that 

of energy density at the edge of the system with respect to 

the entropy density is s = ce'^/'*, we can express the en- 

the energy density at the center. In this way, the function 

tropy change over rapidity from the energy density for- 

/{Rtp) must be less than unity for any angle </>, as the en- 



mula ( 11 ) as: 

ergy density at the center is higher than at the edge. This 

modification from the original Landau assumption plays ~ '^^o ^^P[^iy+ + J/- ^ \/j/+J/-)] ^^gj^^ • ^^^^ 
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aTA{r 



Inserting the solution for the FO time equation (17 1 into 



1-1 



A 



the entropy equation above and assuming that the num- 
ber of produced particles is directly proportional to the Here and everywhere else the vector r starts at the center 
entropy, dN cx dS, one can obtain the number of particles of the almond shaped system of interest. The thickness 
for different rapidities at a fixed angle (j). However, the functions arc then expressed as: 
function f{R,p) and initial thickness A is still needed. 



4 Initial conditions and results 

As is natural for hydrodynamics, the initial state is based 
on the predictions from other models. In this case the 
widely accepted and analytically simple Wounded Nucleon 
(WN) model [T9] will be used to parametrize initial con- 
ditions. It is based on the Woods-Saxon nuclear density 
parametrization ^20^, as follows: 



TAir,q^) =TA{x-b/2,y) ^ J dz p^(r), 
Tsir, 0) = Tb{x + 6/2, y) ^ J dz psW, 



using a Woods-Saxon parametrization ( 19 1 with Ra = 
1.12^1/3-0.86^-1/3 [fm], d = 0.54 [fm] andno = 0-17 fm~^. 

Now assuming that the energy density is proportional 
to the WN density ,21j e(r, (p; b) cx nwN{r, <p; b), the func- 
tion f{R(jj) can be obtained. It is, by definition, the ratio 
of energy density at the edge of the system to the energy 
density at the center, for a fixed impact parameter b and 



l + exp(E^' ^^^^ i^eadsas: 

which is continuous and can be straightforwardly con- , nvKAr(i?0, 6) — min(nvKJv(-R(;i) 0; 6)) , , 

^^""^^^ nw.(0,0;6) ■ (2°) 

nected to the Landau equations. The main requirement for 

the initial conditions and new function f[R^) is that for ^he radius of the system, R^, is dependent on the an- 

i 1 11- • L o j-i li J- u 1 Kle (b and is obtained from the geometry of two overlap- 

the central collision case, o = 0, the result must be equal *= ^ t, j i- 

to the original Landau one. Moreover, the WN model con- Pi^^S ^^^^^^^^ Rl + R^bcos4> + - R\ ^ The term 

nects the impact parameter b with the number of partic- ^HnwN^Rc^, 0; h)) is a minimal density at the edge of 

, ,r J ,1 u r T.- 1 the system and is used in order to have the original Lan- 

ipatmg nucleons, impart, and the number ot binary coi- ^ 

V • T\T 1 • • -4^1, 4^1, -4-1 dau solution for i?A ~ Ra, so that URa) — 0. Because 

lisoms, Ncoii, making comparisons with the experimental ^ ' j \ 

data easy Landau originally used a "sharp sphere" picture, assuming 

The density of wounded nucleons in the transverse ^^^^ ^'^^'"Sy density is zero at the edge of the system 



plane and in polar coordinates, (r, 0), can be obtained ^^^^^ ^^^^^^S lllf' ^^^^ Woods-Saxon model 



by: 



the density ( 19 ) at the edge of the nucleus at r = Ra is 



1 - 



B 



not zero, so the minimal value is subtracted. Now the ac- 



celeration term in ( 14 ) is the same in central collision and 



y'Vbeam 
-3 -2 
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in Landau, but for peripheral collisions the acceleration 
does depend on the angle </>. 

Finally, to calculate the elliptic flow ([2| one should 



0.07 
0.06 
0.05 



merge equations (10 17 18 201 and the initial longitu- 
dinal thickness A, which for the peripheral collisions is 
expressed as: 



Z\(0) = k^Ra/j, 



(21) 
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Cu-Cu -i/s=200 QeV. b=3fm 
Cu-Cu Vs=62.4 QeV, b=3fm 
Au-Au Vs=200 GeV, b=6fm 
Au-Au Vs=62.4 GeV, b=6fm 



Au-Au Vs=200 GeV 10-40% STAR h 
Cu-Cu Vs=200 GeV 0-40% PHOBOS i 
Cu-Cu Vs=62.4 GeV 0-40%. PHOBOS i 
Au-Au Vs=200 GeV 0-40% PHOBOS i 



4 



where — \/nwN{R<in b) / ra.&-K{nw n{R(P, <P\i'))- It nicans 

that the longitudinal expansion ([5| starts with an az- 

imuthally asymmetric initial thickness, which is wider where 

the initial nuclear density is higher. The term maxfnvi/ AT (-R(A, ,, n , , / , n , ^ a 

Fig. 1. Elliptic now dependance on (pseudo) rapidity tor An 

is used in order to have K^{b — 0) — 1 for the cen 



-3 -2 
^'Vbeam 



tral collision case. The energy dependance on the ini- 
tial state is found via the Lorentz gamma factor, 7 = 
y/sN~N /2'mproton- Thc obtained results are shown in fig- 
ure ([1]) for reactions of Au-Au at 6 = 6fm and Cu-Cu at 



An and Cu-Cu for two different collision energies ,/sJrN = 62.4 
GeV and 200 GeV and for different centralities as indicated on 
the figure. The model outcome is on the top, while data from 
RHIC experiments [11] and [12] are on the bottom. 



b — 3fm for two different beam energies, as observed at 

T^TTT^ ^ -1 ;i ; .1 111 1- , ,1 ouW for partlclc multiplicity spectra, but for elliptic flow 

RHIC. One can easily see that the model does predict the i- r- j t- i 

1 ,. 1 1 ^ 1 , 1 11 i ii and elliptic flow longitudinal scaling as well, 

longitudinal scaling, but does not correspond well to the t- a a 



experimental data. The latter can be explained by the lack 
of realistic Freeze-Out dynamics and hadronic rcinterac- 
tions in thc model. Moreover, the initial state is chosen to 
be as simple and transparent as possible, in order to show 
pure hydrodynamic effects for a system with transverse 
asymmetry and does not include initial state fluctuations. 
The initial energy density distribution /(-R^) and initial 



The generalized solution of Landau hydrodynamics in- 
corporates transverse asymmetries into the initial config- 
uration and can be compared to a wide amount of experi- 
mental data, giving deeper insight into the early dynamics 
of the system. 

Acknowledgement K.T. acknowledges the support 



longitudinal thickness might be based on more so- from the Baltic- American Freedom Foundation and Van- 

phisticated models, but for now we can conclude that Lan- derbilt University Relativistic Heavy Ion Group for valu- 
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